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The mean rate of energy dissipation (e) in turbulence is often written in the form of (e) = 
C'ti(u^)'^/^/i„, where the root-mean-square velocity fluctuation (u^)^/^ and the velocity correla- 
tion length Lu are parameters of the energy-containing large scales. However, the dimensionless 
coefficient C„ is known to depend on the flow configuration that is to induce the turbulence. We de- 
fine the correlation length L^2 of the local energy w^, study C„2 — {e}L^2 / (u^)^^^ with experimental 
data of several flows, and find a possibility that C„2 does not depend on the flow conflguration. Not 
Lu but i„2 could serve as the typical size of the energy-containing eddies, so that (m^)^/^/Lu2 is 
proportional to the rate at which the kinetic energy is removed from those eddies and is eventually 
dissipated into heat. 



I. INTRODUCTION 

Since the kinetic energy of turbulence is transferred from large to small scales and is eventually dissipated 
into heat, the mean rate of energy dissipation per unit mass (e) — j{(dxiVj + dx^Vi)'^) /2 is independent 
of the value of the kinematic viscosity v. The rate is instead determined by parameters of the large scales: 

{e)^cS!^. (1) 

^: 

Here C is a dimensionless coefficient, (u^y^'^ is the root-mean-square fluctuation of the velocity component 
Vi in some direction Xi [i = 1, 2, or 3), and the length L represents the large scales or equivalently the sizes 
of the energy-containing eddies. The energy of such eddies is of the order of {v?)- Their time scale is of the 
order of L/ {u^)^/"^. As a result, is of the order of the rate at which their energy is transferred to 

the smaller eddies. This rate is in turn equal to (e). 

Historically, Eq. (1) was found by Taylor^ in 1935. Dryden^ in 1943 and Batchelor'^ in 1953 found that 
Eq. (1) is not inconsistent with experimental data of grid turbulence. Thereafter, Eq. (1) has been used in 
various studies of turbulence, some of which were already described by Landau and Lifshitz* in 1959. 

The length L is traditionally defined as the correlation length L„ of the turbulent velocity u,"^'^ which is 
based on the two-point correlation (u(x -I- r)u{x)): 

T _ {u{x + r)u{x))dr 

However, by using experimental and numerical data of various flows, Sreenivasan^^'' found that Cu = 
depends on the flow conflguration. That is, C„ is determined by the initial condition, boundary 
condition, external force, and so on of the turbulence. For the same flow conflguration, C„ was almost a 
constant of order unity if the Reynolds number was high enough to separate the energy-containing large 
scales from the energy-dissipating small scales. These flndings have been conflrmed experimentally^^^^ and 
numerically. ^^^^^ 

The dependence of C„ on the flow conflguration implies that (u^)'^/^/L„ is not exactly proportional to the 
rate at which the energy is removed from the energy-containing eddies. Although (m^) is not proportional to 
the total energy J^ii'^^f) turbulence is not isotropic, this is not the sole reason because C„ is different 

even among isotropic flows. ^'''■^"^"^ It is hence concluded that L„ is not exactly proportional to the typical 
size L of the energy-containing eddies. 
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TABLE I. Experimental conditions and turbulence parameters of grid turbulence G1-G5, boundary layer B1-B6, and 
jet J1-J6: incoming flow velocity U^t, measurement position a;„t and Zwt, sampling frequency /, moan stroamwise 
velocity U, kinematic viscosity v, mean rate of energy dissipation (e) = 15:y{{dxv)^) /2, root-mean-square velocity 
fluctuations (w^)^''^ and (v^)^''^, correlation lengths Lu, L^, L^2, and L„2 in Eq. (2), and Reynolds numbers ReA„ = 
{uy/^Xu/u and Rex„ = {vY/^X,/iy for A„ = [2{u^)/{{d^v)y^/^ and A. = l{v^}/(id:,v)y^/^. 
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This conclusion is important because Lu has been used as a representative of the large scales, not only in 
the study of Eq. (1) but also in many other studies. It is desirable to explore a more suitable definition of 
L. Pearson et al.^ defined L as a scale that corresponds to the peak of the so-called prcmultiplicd energy 
spectrum, but their definition has turned out to retain the dependence of C on the fiow configuration.^" We 
instead define L as the correlation length L„2 of the local energy u^: 

The expectation is that the fluctuations of represented by L„2 could be related to the energy-containing 
eddies, which arc dcflned to possess the mean energy (m^). By using experimental data of several flows of 
fully developed turbulence, C„2 = (£)L„2/(w^)'^/^ is studied over a range of the Reynolds number. 



II. EXPERIMENTAL DATA 



The experimental data used here are those of the longitudinal velocity u and the lateral velocity v in grid 
turbulence (G1-G5), boundary layer (B1-B6), and jet (J1-J6), among each of which the flow configuration 
was the same but the Reynolds number spans some range. While Bl, B4, B6, Jl, and J6 are used here 
for the first time, the others were used in our past works. ^^'^^ The experimental conditions and turbulence 
parameters are summarized in Table I. 

Especially in an estimation of {e), uncertainties are unavoidable. ^'^'^ To avoid any of the resulting bias, 
the data were obtained and processed in the same manner. The details are described below. For the reader 
not interested in them, we suggest to skip to Sec. III. 



A. Experiments 

The experiments were conducted in a wind tunnel of the Meteorological Research Institute. We adopt 
coordinates a;wt, 2/wti and Zwt in the streamwise, spanwise, and fioor- normal directions. The corresponding 
flow velocities are U + u, v, and w. Here U is the average while u, v, and w are the fluctuations. The origin 
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r/rj r/r] 

FIG. 1. Two-point correlations of u, , v, and as a function of scale r in grid turbulence G3 (circles), boundary 
layer B3 (triangles), and jet J3 (squares). While the scale r is normalized by the Kolmogorov length rj, the correlations 
are normalized by their values at r = 0. The arrows indicate L^, L^2, L„, or L^2. 

Xy/t = ^wt = -^wt = m is on the floor center at the upstream end of the test section of the tunnel. Its size 
was Sxv,t = 18 m, Sy^jt = 3 m, and Sz^jt = 2 m. The cross section Sy^jt x <5^wt was the same upstream to 
a^wt = -4 m. 

To measure U + u and v, we used a hot-wire anemometer. The anemometer was composed of a constant 
temperature system and a crossed-wire probe. The wires were made of platinum-plated tungsten, 5 /xm in 
diameter, 1.25 mm in sensing length, 1mm in separation, oriented at ±45° to the strcamwise direction, and 
280 °C in temperature. Although the spatial resolution of our probe was not high, this is not serious as 
demonstrated in Appendix. 

For the grid turbulence G1-G5, we set a grid at x„t = — 2m across the flow passage to the test section 
of the wind tunnel. The grid had two layers of square rods, with axes in the two layers at right angles. The 
cross section of the rod was 0.04 x 0.04 m^. The spacing of the rod axes was 0.20 m. We set the velocity 
of the incoming flow to be C/wt = 4-20 ms~^. The measurement was on the tunnel axis, ywt = Om and 
z^t = 1.00 m. 

For the boundary layer B1-B6, roughness blocks were set over the entire floor of the test section. The 
block size was 5xy^ = 0.06 m, Jj/^t = 0.21m, and 5zy,t = 0.11m. The spacing of the block centers was 
Sx^t = %wt = 0.50 m. We set the velocity of the incoming flow to be J/wt = 2-20 ms~^. The measurement 
was in the log- law sublayer at .x„t = +12.5 m and j/„t = m, where the boundary layer had the displacement 
thickness of 0.2 m and the 99% velocity thickness of 0.8 m. 

For the jet J1-J6, we set a contraction nozzle. Its exit at Xwt = — 2m was rectangular with the size of 
5ywt = 2.1m and Sz„t = 1.4 m. The center was on the tunnel axis. We set the flow velocity at the nozzle 
exit to be C/wt = 4-41 ms^^. The measurement was at a;wt ~ +15.5 m and y^t = Om. 

The signal of the anemometer was linearized, low-pass filtered, and then digitally sampled. We determined 
the sampling frequency / as high as possible, on the condition that high-frequency noise was not significant 
in the energy spectrum. The filter cutoff was at f/2. We obtained a long record of 1 x 10* data for / < 20 kHz 
or of 4 X 10* data for / > 20 kHz. Since the variation of air temperature was ±1 °C at most, the kinematic 
viscosity v was assumed to have been constant. 

B. Data processing 

The temporal fluctuations u{t„t) and v{t„t) along time t^t were converted into spatial fluctuations of 
the longitudinal velocity u{x) and the lateral velocity v{x) along position x by using Taylor's hypothesis, 
X = —Uty,t. Hereafter, the average (•) is taken over the position x. 
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FIG. 2. Dimensionless coefficients C„ = {e)Lu/{u^)^^^ and C„2 = {e)L^i / {u^)^^'^ as a function of ReA„ and Cv = 
{e}Lv/ (v^)^^^ and C„2 = {e)Ly2 / {v'^)^^^ as a function of ReA„ in grid turbulence G1-G5 (circles), boundary layer 
B1-B6 (triangles), and jet J1-J6 (squares). The solid curve is a fit for C„2 oc Re^^''^ or for C^2 oc Re^y^. 



The mean rate of energy dissipation (e) was calculated as lf>u{{dxvY) /2 instead of usual l'ov{{dxuY) 
by assuming local isotropy, {{dxv)'^) = 2{{dxu)'^). This is because the measurement of the v velocity 
was more reliable. The two wires of the crossed-wire probe individually respond to all the m, w, and w 
velocities. Since the measured u velocity corresponds to the sum of the responses of the two wires, it is 
contaminated with the w velocity especially at smallest scales. Since the measured v velocity corresponds 
to the difference of the responses, it is not contaminated with the w velocity. We have ensured the local 
isotropy by comparing u{x + r) — u{x) with v{x + r) — v{x) along small r. The derivative dxV was estimated 
as [Sv{x + 5x) — 8v{x — 6x) — v{x + 26x) + v{x — 26x)]/126x with the sampling interval Sx = U/f. We also 
calculated the Reynolds number ReA for the Taylor microscale A defined as A„ = {{u^) / {{dxu)'^)Y/^ or as 
A« = [{v^) / {(dx^y)]^^^ , where the local isotropy was again assumed to estimate {{dxu)^} from {{dxv)'^)- 

The correlation lengths in Eq. (2) were calculated for both the u and v velocities. Usually, the two-point 
correlation is integrated only up to the scale of its first zero crossing. This limit is intended to avoid 
statistical uncertainties at the larger scales, albeit at the expense of any information there. Since our data 
records are long, we integrated the two-point correlation beyond the scale of the first zero crossing. The 
convergence of the integration has been ascertained by changing its limit. Figures 1(a) and 1(b) show the 
two-point correlations {u{x + r)u{x)) and {[u'^{x + r) — {u^)][u^(x) — (u^)]). The correlation decays faster 
than the u correlation. As a result, £^2 is smaller than L„ (arrows). This is also the case between the v and 
correlations in Figs. 1(c) and 1(d), although i„2 is larger than L„ in the boundary layer and jet where 
the V correlation is negative at large r. For the same flow configuration, each of L„, Ly, Ly2, and Ly2 has 
almost the same value. 



III. RESULTS 

Figure 2 shows Cu = {e)Lu/ {v^)'^^^ and C„2 = (£)L„2/(u^)^/^ as a function of ReA„ and C„ = 
{e)Ly/ {v"^}^/^ and Cy2 = (e)L^2/(?;^)'^/^ as a function of ReA„. While the sequences of C„ and Cj, do 
not align among the grid turbulence, boundary layer, and jet in Figs. 2(a) and 2(c), the sequences of C„2 
and Cy2 do align in Figs. 2(b) and 2(d). There is a possibility that C„2 and Cy2 are at least approximately 
independent of the flow configuration. As the typical size L of the energy-containing eddies, Z/„2 and Z/^2 
are preferable to L„ and L-^. 

The sequences of C„2 in Fig. 2(b) appear to align better than those of Cy2 in Fig. 2(d). It might follow 
that i„2 is preferable to Ly2 , but any conclusion awaits studies of turbulence that is much more anisotropic. 
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FIG. 3. (color) Normalized variances {{u% - {u\)Y) / {{u^ - {u^)Y) as a function of R/L^2 and {{v% ~ (v\)Y) / {{v'^ - 
(v^))^) as a function of R/L^i in grid turbulence G1-G5 (circles), boundary layer B1-B6 (triangles), and jet J1-J6 
(squares). The solid line indicates the relation of Eq. (4b). 



Since ReA is on a log scale, it is emphasized in Fig. 2 that Cu , C'„2 , Cy , and C„2 are not yet constant but 
continue to decrease with an increase in Re^. The same trend is observed in results of the past works. ^'^"^^'^ 
Such values of Rca are not high enough for complete separation of the large scales from the small scales, 
so that even the energy-containing eddies undergo some energy dissipation in addition to the energy transfer 
oc to the smaller eddies.^'' For convenience to the study in Sec. V, we approximate our results as 

oc Re^" and C^2 cx Re^" with a — 1/2 (solid curves). The exponent a is nevertheless not constant and 
should become negligible at the higher value of the Reynolds number ReA.^""'^^ 



IV. DISCUSSION 



Having found a possibility that L„2 could serve as the typical size L of the energy-containing eddies, its 
reasoning is discussed here. The data record of u{x) is divided into segments with length R. For each of the 
segments, the center of which is tentatively defined as Xc, the energy v? is coarse-grained as 



1 



R 



+R/2 



'^r{xc) ^ S / {xc + x)dx. (3) 



i?,/2 



If u{x) is homogeneous, the mean square of around its average (it|,) 



2 



R 



((4 - {nl) f) = J^J^ (R- r){[n'i^ + r) - {u')][uHx) {u')])dr, (4a) 

where (•) denotes both the averages over the positions and over the segments. Also if the correlation is 
negligible at r 3> i„2, Eqs. (2b) and (4a) yield a relation of {{uj^ — (w^))^) to the correlation length 

((4 - (4))') = ^(("' at i?»L„2. (4b) 

We relate {{uj^ — (w^))^) to the typical size L of the energy-containing eddies that have been defined to 
possess the mean energy (u^). The segment with length R is divided into subsegments with length nL <C R- 

^ R/nL 

^R = ^ "«L(2;m), (5a) 

m— 1 

where Xm is the center of the mth subsegment. If n = 1, the subsegments are just the energy-containing 
eddies with the mean energy {uj^) = {u^). The adjacent eddies might be correlated with one another, but 
such a correlation is negligible if we collect a sufficient number of them, n — N ^ 1. Since the mean square 
of Y^fyl=i' ''^NLi^rn) arouud its average is {R / N L) {{uj^ — we obtain from Eq. (5a) as^^ 

((4-(4))') = ^(«l-«l))') at i?»7VL»i. (5b) 



6 



Then, Eqs. (4b) and (5b) yield 



L.=,L for ^ = i^((!fk^(!fk))!) (5c) 



If some fixed values of A'' and 7 are used to determine i, it is proportional to L^i . The values of A'' and 7 are 
not determined without a further assumption, -'^^ but here wc discuss them a little more. Figure 3(a) shows 
((^H ~ (^fl))^) our experiments. They agree with Eq. (4b) at i? > 10^L„2 (solid line). Since Eq. (4b) 
corresponds to 7 = 1 in Eq. (5c) if we set R = NL, N ~ 10^ is enough for 7 ~ 10'^ to have L„2 = 7L 
irrespective of the flow. The result also ensures that our estimates of L„2 are reliable. As observed in 
Fig. 3(b), the same discussion applies to L„2. 

The above discussion does not apply to Z/„ or Ly. It is true that L„ is related to the variance of the 
coarse-grained velocity Uii{xc) = u{xc + x)dx/R in a manner similar to Eq. (4b): 
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{Ul) = ^{u') at (6) 

However, such a fluctuation of Ur is not related to the typical size L of the energy-containing eddies so 
far as they are defined to possess the mean energy (u^). The reason is {u^) > (f/^) at any R. We know 
of no theory that prefers i„ or Ly as the typical size L. They might rather suffer from too large scales 
(Fig. 1). For isotropic turbulence, Batchelor^ pointed out that Ly, oc J^k~^E{k)dk does not represent 
the part of the three-dimensional energy spectrum E{k) that makes the major contribution to the energy 

OC E{k)dk. 



V. CONCLUDING REMARKS 



For the mean rate of energy dissipation written in the form of {e) = C {v?)'^^'^ / L, it is traditional to define 
L as the velocity correlation length L„. However, Cy = {s)Ly/ (u^)^^^ depends on the flow configuration 
that is to induce the turbulence. We have defined L as the correlation length L„2 of the local energy u^, 
studied C„2 = {e)Ly2 /{u^}"^/^ for several flows, and found a possibility that C„2 does not depend on the flow 
configuration in Fig. 2. Not L„ but L„2 could be proportional to the typical size L of the energy-containing 
eddies, so that (u^)^/^/L„2 is proportional to the rate at which the kinetic energy is removed from those 
eddies and is eventually dissipated into heat. The energy-containing eddies, if defined to possess the mean 
energy (m^), correspond to spatial fluctuations of represented by Ly2 rather than to those of u represented 
by Ly. This is also the case for the correlation lengths Ly and Ly2 of the velocity v and of the local energy 
v\ 

The lengths Ly2 and L„2 describe the variances of u% and at large R (Fig. 3). In fact, our past 
work^*' has demonstrated that L„2 and Ly2 describe their entire distributions. We consider Ly2 and Ly2 as 
important lengths of the energy-containing scales. 

We have an implication from the K-e model^° of the numerical simulation, which calculates the eddy 
viscosity at each position as i't = cK"^ / {e). Here = ^j{vf ) /2 and (e) are ensemble averages. Defined as 
c=cC for VT = cKL and C = {e)L/K^^'^, the constant c = 0.09 applies to various flows. Then, a constant 
value of c could also apply to the various flows if L could be defined as the correlation length of some local 
energy so that C is independent of the flow configuration. 

The existing discussions often assume that C = {e)L/ {u"^)^^^ is independent of the Reynolds number 
Re^.^'* They have to be corrected at Re^ < 10^, where we observe C oc Re^" (Fig- 2). The large- 
scale Reynolds number {u^y^^L/u is corrected as oc Re^~". The number of degrees of freedom per unit 

volume {L/riY, where 77 = {v'^ /{e)^/^ is the Kolmogorov length, is corrected as oc Re^^^~^". Also if, say, 
Loitsyansky's invariant holds as oc (u^)L^ in decaying isotropic turbulence, dt{v?) oc —(e) leads to the decay 
law (m2) oci-(i"-5")/(7-5"). 

Since our data set is not extensive, we have only proposed a possibility that L„2 or Ly2 could serve as L. 
This has to be examined in future with a more extensive set of experimental or numerical data especially at 

high Reynolds numbers, Rca > 10'^. Although we are not so confident that L„2 or i„2 does serve as i, it 
is important to persistently explore a suitable definition of L or of the typical size of the energy-containing 
eddies. 



300 1000 



300 1 000 



FIG. 4. Dimensionless coefficients Cu = {e)Lu/ {v?)^^^ and C„2 = (e)L„2/(M^)^/^ as a function of ReA„ in grid 
turbulence G1-G4 (circles), boundary layer B2-B5 (triangles), and jet J2-J5 (squares). The filled symbols denote 
data obtained with a single-wire probe. The open symbols denote data obtained with a crossed-wire probe. 
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Appendix: SUPPLEMENTARY EXPERIMENTS AND THEIR RESULTS 

The spatial resolution of our crossed-wirc probe was not high. We used a single-wire probe to conduct 
supplementary measurements of the streamwisc velocity u for the configurations of G1-G4, B2-B5, and J2- 
J5. The wire parameters were the same as those of the crossed-wire probe. We estimated (e) as lf>v{{dxu)'^). 
It has been ascertained that the values of —{{dxu)^)/{{dxu)'^)^^'^ agree with those of experimental and 
numerical data in the literature, 0.5-0.6 for Hex ~ 10^-10^. The values smaller by 0.1-0.3 were obtained 
with the crossed-wire probe. 

Figure 4 compares the results for the single-wire probe (filled symbols) with those for the crossed-wire 
probe (open symbols). They agree with each other. The exception is C„ — {e)Lu/ (u^)^^"^ in the jet where Lu 
is different by 15-25% (squares). This is because the u velocities measured with the single- and crossed-wire 
probes were not identical. To the traverse system of the wind tunnel, the probe was set in a different manner. 
The contamination with the v and w velocities was also different (see Sec. II B). Since the alignment of the 
sequences of C„2 = {e)L^2 / {u'^)^/^ is again observed in Fig. 4(b), the spatial resolution of our crossed-wire 
probe is not serious to the estimates of {e), albeit serious to those of (e^), {e^}, and so on, and hence is not 
serious to our conclusion that C„2 is independent of the flow configuration at least as an approximation. 
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